Small time asymptotics in local limit theorems for 
Markov chains converging to diffusions. * 

Valentin KONAKOV f 
2nd February 2008 



We consider triangular arrays of Markov chains that converge weakly to a diffusion 
process. Local limit theorems for transition densities are proved. The observation time 
[0, T] may be fixed or lim n _+ooT = 0, where nh = T and h is a mesh between two 
neighboring observation points. 

1991 MSC: primary 62G07, secondary 60G60 

Keywords and phrases: Markov chains, diffusion processes, transition densities, Edge- 
worth expansions 



*This research was supported by grant 436RUS113/845/0-1 from the Deutsche Forschungsgemein- 
schaft and by grants 04-01-00700 and 05-01-04004 from the Russian Foundation of Fundamental 
Researches. The author worked on the paper during a visit at the Laboratory of Probability Theory 
and Random Models of the University Paris VI in 2005/2006. He is grateful for the hospitality during 
his stay. 

^Laboratory of Probability and Statistics, Central Economics Mathematical Institute, Academy of 
Sciences, Nahimovskii av. 47, 117418 Moscow, Russia, E mail: kv24@mail.ru 



1 



1. Introduction and main results 



Let n > 2 and h > be such that T — nh < 1. Suppose that q (t,x, •) , (t,x) G 
[0, 1] x l d is a given family of densities on M d and m is a function from [0, 1] x ~R d into 
M. d . We shall impose the following conditions 

Al J Rd yq(t,x,y)dy = 0, < t < 1, x G R d . 

A2 There exists positive constants a* and a* such that the covariance matrix a (t, x) = 
f R d yy T Q (t, x, y) dy satisfies 

a* < 9 T a(t,x)9 < a*, 

for all = 1 and t G [0, 1] x G M d . 

A3 There exists a positive integer S" and a real nonnegative function vp (y) , y £ W d 
satisfying sup ygR d ip (y) < oo and f Rd \\y\\ ip (y) dy < oo, with 5 = 2c?S" + 4, such 
that 

\D^q(t,x,y)\ < i e [0,1], |H =0,1,2,3,4 

and 

\r%q(t,x,y)\<ij>(y), te [0,1], H =0,1,2. 

Furthemore, for all x G M d it holds 




q(t,x,y) - q(t',x,y)\dy -> 0, 



as |t-f| ->0. 

Bl The functions m(t,x) and <r(t,x) and their first and second derivatives w.r.t. t 
and x are continuous and bounded uniformly in t and x. All these functions are 
Lipschitz continuous with respect to x with a Lipschitz constant that does not 
depend on t. Furthemore, D v a(t 1 x) exists for \v\ = 6 and is Holder continuous 
w.r.t. x with positive exponent and a constant that does not depend on t. 

Consider a family of Markov processes in M. d of the following form 

(1) X k+lth = X kth + m (kh, X kth ) h + y/h£ k +i,h, x o,h = x G R d , k = 0, n - 1, 

where (£i,/i) i=1 n is an innovation sequence satisfying the Markov assumption: the con- 
ditional distribution of ^k+\,h given — x^, X 0y h = xq depends only on = x k 
and has conditional density q (kh, x k , •) . The conditional covariance matrix correspond- 
ing to this density is a(kh, x k ). The transition densities of (X it h) i=1 are denoted by 
p h (0,kh,x,-). 
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We shall consider the process (JTJ) as an approximation to the following stochastic 
differential equation in R d : 



dY s = m (s, Y s ) ds + A (s, Y s ) dW s , Y = xeR d , s G [0, T], 
where (W s ) s>0 is the standard Wiener process and A is a symmetic positive definite 



d x d matrix such that A (s, y) A (s, y) T = a (s, y) . The conditional density of Y t , given 
Y = x is denoted by p (0, t, x, •) . 

Konakov and Mammen (2000) obtained a nonuniform rate of convergence for the 
difference pt (0, T, x, •) — p (0, T, x, •) as n — > oo in the case T x 1. It is the goal of the 
present note to get an analogous result in the case T — o(l) . The following theorem 
is our main result. 



Theorem 1 Let h > 0, n > 2, T = nh. Assume (A1-A3) and (Bl). Then, as n — ► oo 
anc? T — > 0, 



sup 



<2v/t - x ) Ipa (°> T > x ,y)-p (°> x, y)| = O (n 1/2 ) , 



where the constant in O (•) does not depend on h and 

2S'-2 



Qs (u) = 5 d 1 + 



2. Parametrix method for difFusions 

For any s 6 (0, T), x, y G M d we consider an additional family of "frozen" diffusion 
processes 

= m (t, y) dt + A (t, y) dW t , F s = x, s<t<T. 

Let (s,f, x, •) be the conditional density of Y t , given y s = x. In the sequel for any z 
we shall denote p(s, t, x, z) = p z (s,t,x,z) , where the variable z acts here twice: as 
the arument of the density and as defining quantity of the process Y t . 

The transition densities p can be computed explicitly 



p(s,t,x,y) = (2tt) d/2 (deter (s,t,y)) 1/2 



x exp (y - x - m (s, t, y)) T a 1 (s,t,y) (y - x - m (s, t, y)) 

where t t 

o~(s,t,y)= / a(u,y)du, m(s,t,y)= / m(u,y)du. 
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We need the following kernel 



1 d 

H(s,t,x,y) = -^2(aij(s,x) -aij(s,y)) 



d 2 p(s,t,x,y) 

+ ^2 ( m * ( s ' x ) - m i ( s ' y)) — ^r 2 — ■ 

1=1 

Introduce the convolution type binary operation ® : 

/ ® g(s,t,x,y) = du f (s, u, x, z) g (u, t, z, y) dz. 

J s JR d 

fc-fold convolution of H is denoted by H^ k > . The following results are taken from Kon- 
akov and Mammen (2000). 

Lemma 2 Let < s < t < T. It holds 

oo 

p(s,t,x,y) = ^2p(g> H (r) (s,t,x,y) . 

r=0 

Lemma 3 Let < s < t < T. There exist constants C and C\ such that 

\H(s,t,x,y)\ < Cip~ 1 (j) C , p (y - x) 

and r 

\p® (s, t, x, y) | < C{ +1 r * ] - fo, p (y-x), 

where p 2 = t — s, (pc, P (u) = P d( Pc ( u /p) an d 

exp (-C ||m|| 2 ) 



<Pc (u) 



J exp (-C \\v\\ ) dv 



3. Parametrix method for Markov chains 

For any < jh < T, x, y G M. d we consider an additional family of "frozen" Markov 
chains defined for jh < ih < T as 

(2) X i+ i th = X iA + m (ih, y)h + y/h£ i+llh , X jih = x eR d , j < i < n, 

where £j+i t h, ■ £n,h is an innovation sequence such that the conditional density of £i+i,h 
given Xi t h = Xi,...,X ^ h = x equals to q(ih,y,-). Let us introduce the infinitesimal 
operators corresponding to Markov chains JTJ) and Q respectively, 

L h f(jh,hh,x,y) 

= h~ X ( / p h (jh, (j + 1) h, x, z) f ((j + 1) h, kh, z, y) dz - f ((j + 1) h, kh, z, y) 



and 



L h f (jh,kh,x,y) 

= h- 1 (Jp y h (jh, (j + 1) h, x, z) f ((j + 1) h, kh, z, y)dz-f ((j + 1) h, kh, z, y)J , 

where p\ (jh,j'h, x, •) denotes the conditional density of Xj> th given Xj h = x. As before 
for any z denote pih(jh, j'h,x, z) = p z h (jh, j'h,x, z) , where the variable z acts here 
twise: as the arument of the density and as defining quantity of the process X^h- For 
technical convenience the terms / ((j + 1) h,kh, z,y) on the right hand side of Lhf 
and Lhf appear instead of / (jh, kh, z, y) . 

In analogy with the definition of H we put, for k > j, 

H h (jh, kh, x, y) = yL h - L h ^j p h (jh, kh, x, y) . 
We also shall use the convolution type binary operation ®h '■ 

k-l . 

9®hf (jh, kh,x,y) = y^h g (jh, ih, x, z) f (ih, kh, z, y) dz, 

where < j < k < n. Write g ® h H ( h 0) = g and g ® h = (g ®h H^' 1 ^ ^ H h 

for r = 1, ...,n. For the higher order convolutions we use the convention Yl\=j = f° r 
/ < j. One can show the following analog of the "parametrix" expansion for ph [see 
Konakov and Mammen (2000)]. 

Lemma 4 Let < jh < kh < T. It holds 

k-j 

p h (jh, kh, x, y) = ^Ph ®h H { h r) (jh, kh, x, y) , 

where 

Ph (jh, jh, x, y) = p h (kh, kh, x,y) = 5 (y - x) 
and 8 is the Dirac delta symbol. 



4. Auxiliary statements 

In this section we collect several bounds to be used later in the proofs. Throughout 
the proofs C,Ci,C 2 , ... denote generic constants possibly different in different places. 

We give some tools that are useful for comparison of the expansions of p and ph 
(see Lemmas [2} and 0] respectively). Since p and ph are written in terms of pft and p^ 
it is essential first to bound the difference p^ — p y ■ For this we use some bounds from 
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Konakov and Mammen (2000) which are proved for T = 1, but remain true also for 
T — > 0. The following lemmas are slight modifications of Lemmas 3.8, 3.9, 3.11, 3.12, 
3.13 from Konakov and Mammen (2000) and therefore the proofs will be not detailes 
here. 



To simplify the notation let p = ^/ (k — j) h. Denote ( p (z) = p p ( (z/p) , where 

(i+Mr)"' 

CM- ■ 



/(l+IMI 5 - 1 ) -1 



dz 



Lemma 5 Assume Conditions (A1-A3) and (Bl). Then for < j < k < n and all 
x, y e R d it holds 

\p h (jh, kh, x, y) - p(jh, kh, x,y)\ < Ch l/2 p' l C P (y - x) . 



For j = 0,...,k-2, let 

K h (jh, kh, x, y)= (l- l)j p h (jh, kh, x, y) 

and 

M h {jh,kh,x,y) = Sh^YsYl f dd f d6D% q (jh,y,e)(x-yy 

x ^D v x p h ((j + 1) h, kh, x + bQh>l\ y) (1 - 5f . 

If j = k — 1 define 

K h ((k - l)h, kh,x,y) = and M h ((k — 1) h, kh, x, y) — 0. 



Lemma 6 Assume Conditions (A1-A3) and (Bl). Then for < j < k < n and all 
x, y G R d it holds 



\H h (jh, kh, x, y) - K h (jh, kh, x, y) - M h (jh, kh, x,y)\ < Ch 1/2 p 1 ( p (y 
where ( p is defined above. 

Denote £ p (z) = p~ p £ (z/p) , where 



x) 



fU+wz'ir-'Y'dz'' 
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Lemma 7 Assume Conditions (A1-A3) and (Bl). Then for r = 1,2,..., < j < k < 
n and all x, y e M. d it holds 

C r+l p r 



p h ® h Hj[> (jh, kh, x, y) 



< 



r(i + |) 



i P (? - y) 



Lemma 8 Assume Conditions (A1-A3) and (Bl). Then for < j < k < n and all 
x,y <ER d it holds 



k-j 



r=0 



p h {jh,kh,x,y) = \Ph ®h (M h + K h ) {r) ) (jh, kh,x,y) + R 1 (jh,kh,x,y) 



where 



(jh,kh,x,y)\<Ch}' 2 p(i P (y-x) 



Lemma 9 Assume Conditions (A1-A3) and (Bl). Then for < j < k < n and all 
x, y e R d it holds 

k-j 

p h (jh, kh, x, y) = ^2 (p®h (M h + K h ) {r) ^j (jh, kh, x, y) + R 2 (jh, kh, x, y) , 



where 



\R 2 (jh, kh, x,y)\< Ch x l 2 p-% (y-x). 



5. Proof of the main result 

From Lemmas [2] and 01 we get as n — > oo and T = nh — > 



(3) p (0, T, x, y) = (p ® H {r) ) (0, T, x, y) + T~ d ' 2 exp 



r=0 



C\\y-x\ 
T 



(T n e- n ) 



Furthemore, Lemma EH implies that 

n 

(4) p h (0, T, x, y) = {v®h (M h + K h ) {r) ) (0, T, x, y) + £ VT (y-x)0 (n^ 2 ) 

r=0 

Because of (JHJ) and (J3J for the statement of the theorem it remains to show that 

n n 

Y,{P® (0, T, x,y)-J2 (P®h (M h + K h ) {r) ) (0, T, x, y) 

r=0 r=0 

(5) =^ VT (y- x)0 (n- 1 / 2 ). 



For the proof of (JHJ note that 



n n 

(p ® (0, T, x,y)-J2 (v®h (M h + K h ) {r) ) (0, T, x, y) 



(6) 



r=0 

< Si + 1S2 + S3, 



r=0 



where 



Si 



So 



J2(P® H^) (0, T, x, y) - J2 {P®h H^) (0, T, x, y) 

r=0 r=0 

n n 

(P®h H {r) ) (0, T, x,y)-J2 (P®h (M h + ) W ) (0, T, x, y) 

r=0 r=0 
n n 

{v®h (M h + #) (r) ) (0, T, x,y)-J2 (P®h (M h + K h f^ (0, T, x, y) 



r=0 r=0 

For 5i, S2, S3 we will show the following estimates 

(7) S k = Q VT (y-x)0 (h 1 ' 2 ) , k = 1,2,3. 

We shall prove (JJJ) for k = 1. The term Si corresponds to the passage from the con- 
tinuous time to the lattice time. The fact that T — > implies that integrands involved 
in the convolutions (g) and 0^ become asymptotically degenerate and therefore more 
accurate estimates are required than those in Konakov and Mammen (2000). We will 
develop here the details for these bounds. 

We start from the recurrence relations for r = 1, 2, 3, ... 

(p®ff (r) ) (0,jh,x,y)- {p® h H {r) ) (0,jh,x,y) 



(8) 



[(p ® H^) ® H - (p®H( r ~V) ® h H] (0,jh,x,y) 



+ [(p® H^) - (p® h H^)] ® h H (0, jh, x, y) . 



By summing up the identities in (JHJ) from r = 1 to 00 and by using the linearity of the 
operations ® and ®h we get 

(P - /) (0, A x,y) = {p®H -p® h H) (0, j/i, x, y) 



(9) 

where we put 
(10) 



+{p-p d ) ® h H{Q,jh,x,y) 



p d [i 



h,i'h,x,y) = }Xp ®h H (r) )(ih,i'h,x,y). 

r=0 



By iterative application of (jHJ) we obtain 

(P - (0, A x,y) = {p®H -p® h H) (0, j/i, x, y) 
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(11) + (p®H-p® h H)® h <S>{0,jh,x,y), 

where $(z/i, i'h, z, z') = H(ih, i'h, z, z')+H® h H(ih, i'h, z, z')+... = J2T=i H ir) (ih, i'h, z, z') 
By the Taylor expansion we have 

(p®H-p® h H) (0,jh,x,z) 

J -1 f (i+l)h 



V / du I [X{u)-X {ih)} dv 
/ (u — ih)du l X'(ih)dv 

4-n Jih JR d 



i ill r(i+i)A /-i r 

(12) +- V / (u - ih) 2 (1-5) X"(s) | s=Si dvdSdu, 

* i=0 Jih JO JR d 

where A (u) = p(0, u, x, v)H(u,jh, v, z), Sj = Si(u, i, 5, h) = ih + S(u — ih). 
Note that 

X'(ih)dv = / —p(0,s,x,v)\ s=ih H(ih,jh,v,z)dv 

R d J R d OS 

+ [ p(0,ih,x,v)-^-H(s,jh,v,z)\ s=ih dv = [ L f p(0, ih, x, v ) 

JRd OS J Rd 

x(L — L)p(ih, jh,v, z)dv — I p(0,ih, x,v)[(L — L)Lp{ih, jh,v, z) 

Jn d 



—H\{ih,jh,v,z)\dv = / p(0,ih,x,v)Hi(ih,jh,v,z)dv 

Jn d 

(13) +/ p(0,ih,x,v)(L 2 -2LL + L 2 )p(ih,jh,v,z)dv, 
Jn d 

where Hi(s, t, v , z) is defined below in (|2TJ) . We get from (JTJJ) 

y / (u — ih)du I X'(ih)dv = — (p <S>h Hi)(Q, jh, x, z) 
~^ Jih Jn d 2 



(14) +^(p® h A o )(0,jh,x,z), 

where A Q (s, jh,v, z) = [L 2 — ILL + L 2 )p(s, jh,v, z). The direct calculation shows that 

1 d 

A (s,jh,v,z) = - ^ (<? Pq {s,v) - a pq (s,z)){a r i{s,v) - a rl (s,z)) 

p,q,r,l=l 



X 



d 4 p(s,jh,v,z) 
dvpdv q dv r dvi 



+ E 



[a pq (s,v) - a pq (s,z))(m r (s } v) - m r (s,z)) 



p,q,r=l 



(15) 



x- 



d 3 p{s,jh,v,z) 

dVpdv q dv r 



1 d 
+ 5 E 



2 



da r i{s,v)d 3 p{s,jh,v,z) 
<7„ n [S,V) — - — h (< 2), 



(9f p dv q dv r dvi 



where we denote by (< 2) the sum of terms containing the derivatives of p(s,jh, v, z) 
of the order less or equal than 2. Note that for a constant C < oo and any < e < \ 



h 

-(p®fc#i)(0, jh,x,z) 



c,VJh 



[z — X) 



(16) 



-(p® h A )(0,jh,x,z) 



<C{e)h^r {1,2 - e) ct>c^-^ 



First inequality (fT6|) follows from (Bl) and the well know estimates for the diffusion 
density p and for the kernel Hi . The second inequality (fT6|) follows from (Bl), ([T| 
and the following estimate 



§E* 



i=0 



p{0, ih, x, v) — dv 

R d OV q OV r OVi 



§E* 



i=0 



dp(0, ih, x, v) 



R d 



d 2 p(ih,jh, v, z) 
dv r dvi 



dv 



i_1 1 1 

< Ch l - £ V —rz JTh {Z-X)< Ch 1/2 



i=0 



(17) xj-^B(^s)4> CiVJn (z-x). 

Now we shall estimate the second summand in the right hand side of (fT2|) . Clearly 

d 2 d 
= -q^P^i s i x i v ) H ( s dh,v, z) + 2—p(0,s J x J v) 



(18) 



d d 2 



Using forward and backward Kolmogorov equations we get from (jTHI) after long but 
simple calculations 



i J -1 r(i+l)h pi p 



i=0 



ill 



dvdbdu 



( 19 ) / (u-ih) 2 / (1-5) V] / p(0, s,x,v)A k (s, jh,v,z) \ s=Si dvdSdu, 

1 i= o J * Jo k=i jRd 



k=l 
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where 



Ai(s,jh,v,z) = (L 3 — 3L 2 L + 3LL 2 — L 3 )p(s,jh, v, z), 
A 2 = {L 1 H + 2LH 1 ){s,jh,v,z), 

A 3 (s,jh,v,z) = [(L — L)Li + 2(Li - Li)L]p(s, jh, v, z), 



(20) 
and 



(21) 



A 4 (s,jh,v,z) = H 2 {s,jh,v,z). 



Hi(s, t, v, z) = (Li - Li)p(s, t, v, z) 



1 d 

-Y 

2 ^ 



d l aij(s,v) d l a i:j (s,z)\ d 2 p(s,t,v 1 z) 



ds l 



ds l 



dvidvj 



+ d l mi(s,v) _ d l mi(s,z) ^ dp(s,t,v,z) , . ( } 



i=i 



9s' 



9s' 



<9t>, ; 



Using integration by parts and the definition |20l) of A 2 , A 3 and A 4 it is easy to 
that for any < e < 1/2 and for k = 2, 3, 4 



J -1 r(i+l)h pi 

{u-ihf j (1-5)/ p(0,s,x,u)A fc (s, jh,v,z) 



i=0 



dvdSdu 



R d 



(22) ^C^^/^^^-x). 

For fc = 1 we shall prove the following estimate for any < e < | 



J -1 p(i+l)h pi 

/ (u — ih) 2 j (1 — 5) j p(0,s,x,v)Ai(s,jh,v,z)\ s=s ^dvd5du 



i=0 



i/i 



(23) ^(e^r^^c,^ (*-*)■ 

Note that the function Ai(s,jh, v, z) can be written as the following sum 
1 d 

At(s,jh,v,z) = - ^2 ( a ij( s i v ) ~ Vij( s , z ))(vpq(s,v) -apg(s,z))(ai r (s,v) 



-cri r (s,z)) 



i,j,p,q,l,r=l 

d 6 p(s,jh, v, z) 
dvidvjdvpdv q dvidv r ' 4 



d 

1 £ < 

11 



(s,z))(a pq (s,v) 



-a (s z))(m,(sv)-m,(s z)) ^.JM,*) + 3 A , . da^v) 
(24) Ms ,„)- Ms , 2 ))x||£^- + (<4), 

where we denote by (< 4) the sum of terms containing the derivatives of p(s,jh, v, z) 
of the order less or equal than 4. By (Bl) and |{2*1J) it is clear that the estimate for the 
left hand side of (f2*2l) for k = 1 will be the same (up to a constant) as for the following 
sum for fixed p,q,r,l 



( {i+1)H < -M* f\i A\ f fn , &p(s,jh,v,z) 
} (u-ih) / (1-6) p(0,s,x,v) — \ s=Si dvdbdu 

t^Jih Jo Jr* dv p dv q dvidv r 



After integration by parts w.r.t. v p and with the substitution hw — (u — ih) in each 
integral we obtain 

i ill r(i+i)h r 1 r &n 



\ I (u-ih) I (1-6)1 p(0,s,x,v) — \ s=Si dvd6du 

t^Jih Jo Jr* dv p dv q dvidv r 



i=0 

1 

2 



[ {i+l)H ( -h\* f\i ^ f dp(0,s,x,v)d 3 p(s,jh,v,z) 



j=0 



o J Rd dv p dv g dvidv r 

i i-i 



< Ch 3/2 - £ 4> C)VJE (z - x) [ w 2 dw [ (l-6) 1/2 ~ £ d6 [ 

Jo Jo Jo 



8=0 

^ h dt 



Vt[(j-i)h-t] 



1-e 



(25) 



<^j-( 1 / 2 - £ )5(^)0 CiV ^( 



2' -/ro,yjra \^ 

where _B(p, g) is a Beta function and (pc,p {z — x) is defined in Lemma El ■ As we 
mentioned above (23J) follows now from (25J. Bv (TT1. (TT41. (fTT)|l . (gZ| and (TZ3jl we 

obtain for any < e < \ and j = 1, 2, ...n 

(26) \(p®H-p® h H) (0, jfc, *) | < C(e)h 1 / 2 r {1/ ^ £) 4>c,VJh (z ~ x) 

We use now the following estimate for <&(ih, i'h, z, z') proved in Konakov and Mammen 
(2002) 

(27) mih,i'h,z,z')\<C 7 =L= cf> C!s/¥K = rh (z'-z) 

\/vh — in 
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From (jTTf pfijl and (|27jl we obtain 

(28) | (p - p d ) (0, n/i, x,y)| < C^)/^ V^ 2 ^^ (y - x) . 

The last inequality proves for k = l.The terms S2 and S3 can be handled in the 
same way as the terms T 2 and T 3 in Konakov and Mammen (2000). This completes 
the proof of our main result. 

Remark 1. In fact for Si we proved stronger result than the estimate 0. We get 
the following representation 

(p - p d ) (0, T, x, y) = -(p ® h F0 (0, T, x, y) + -(p ® h A,) (0, T, x, 
+ -(p ® h E x ® h $) (0, T, x, y) + -(p ® h A ® h $) (0, T, x, y) 



(29) + J R(0,T,x,y), 
where for any < e < 1/2 

172(0, T, x,y)\ < C(e) (/i 3/2 ~ £ + /m - ^-*) + h^T)<P c VT (y - x) 

= C(e)<f>c,y/r(y- x ) °( h )- 

This representation is useful to obtain a small time Edgeworth type expansions for our 
model (see Konakov and Mammen (2005)) 

Remark 2. If T = const, (without loss of generality we assume T = 1) than we 
can easily avoid the difficulties connected with singularity by splitting the time interval 
[0,1] and by using an integration by parts. For example 

du / p(0,u,x,v) dv= / ...+ / ... 

Jri dVpdVqdvidVr J J 1/2 

du \ p{0,u,x,v) dv 

jRd dVpdVgdVidVr 



(30) + /W ggg$n M)dv 

Jl/2 



R d dvpdvqdvidv r 



and the derivatives in the right hand side of (j30l) are not singular. The representation 
(|2~H| remains true. All summands in the right hand side of (|2~H1) are estimated from 
above in absolute value by Ch(j)c,i (y — x) and for the remainder term R(Q, 1, x, y) the 
following estimate holds 

(31) \R$A,x,y)\<Ch 2 <j> c ,i(y-x). 
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